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Abstract — Quasi-cyclic (QC) low-density parity-check (LDPC) 
codes are an important instance of proto-graph-based LDPC 
codes. In this paper we present upper bounds on the minimum 
Hamming distance of QC LDPC codes and study how these upper 
bounds depend on graph structure parameters (like variable 
degrees, check node degrees, girth) of the Tanner graph and 
of the underlying proto-graph. Moreover, for several classes of 
proto-graphs we present explicit QC LDPC code constructions 
that achieve (or come close) to the respective minimum Hamming 
distance upper bounds. 

Index Terms — Girth, graph cover, low-density parity-check 
code, low-density parity-check matrix, proto-graph, proto-matrix, 
pseudo-codeword, quasi-cyclic code, Tanner graph, weight ma- 
trix. 



I. Introduction 

^""V UASI-CYCLIC (QC) low-density parity-check (LDPC) 
\J codes represent an important class of codes within 
theramily of LDPC codes fTJ. The first graph-based code 
construction that yielded QC codes was presented by Tanner 
in 0; although that code construction was presented in the 
context of repeat-accumulate codes, it was easy to generalize 
the underlying idea to LDPC codes in order to obtain QC 
LDPC codes 0, 0, 13, 0. The simplicity with which 
QC LDPC codes can be described make them attractive for 
implementation and analysis purposes. 

A QC LDPC code of length n = Lr can be described by a 
Jr x Lr (scalar) parity-check matrix that is formed by a J x L 
array of r x r circulant matrices. Clearly, by choosing these 
circulant matrices to be low-density, the parity-check matrix 
will also be low-density. 

With the help of the well-known isomorphism between the 
ring of circulant matrices and the ring of polynomials modulo 
x T — 1, a QC LDPC code can equally well be described by a 
polynomial parity-check matrix of size JxL. In the remainder 
of the paper we will mainly work with the polynomial parity- 
check matrix of a QC LDPC code and not with the (scalar) 
parity-check matrix. Another relevant concept in this paper 
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will be the weight matrix associated with a polynomial parity- 
check matrix; this weight matrix is a J x L integer matrix 
whose entries indicate the number of nonzero coefficients of 
the corresponding polynomial in the polynomial parity-check 
matrix. 

Early papers on QC LDPC codes focused mainly on poly- 
nomial parity-check matrices whose weight matrix contained 
only ones. Such polynomial parity-check matrices are known 
as monomial parity-check matrices because all entries are 
monomials, i.e., polynomials with only one nonzero coeffi- 
cient. For this class of QC LDPC codes it was soon established 
that the minimum Hamming distance is always upper bounded 
by (J + l)! 0, 0, 13. 

In this paper we study polynomial parity-check matrices 
with more general weight matrices by allowing the entries 
of the weight matrix to be 0, 1, 2, or 3 (and sometimes 
larger). This is equivalent to allowing the entries of the 
polynomial parity-check matrix to be the zero polynomial, 
to be a monomial, to be a binomial, or to be a trinomial 
(and sometimes a polynomial with more nonzero coefficients). 
The main theme will be to analyze the minimum Hamming 
distance of such codes, in particular by studying upper bounds 
on the minimum Hamming distance and to see how these 
upper bounds depend on other code parameters like the girth 
of the Tanner graph. We will obtain upper bounds that are 
functions of the polynomial parity-check matrix and upper 
bounds that are functions of the weight matrix. The latter 
results are in general weaker but they give good insights into 
the dependency of the minimum Hamming distance on the 
structure of the weight matrix. For example, for J = 3 we 
show that there are weight matrices that are different from 
the all-one weight matrix (but with the same column and row 
sum) that yield minimum Hamming distance upper bounds 
that are larger than the above-mentioned (J + l)! bound. By 
constructing some codes that achieve this upper bound we are 
able to show that the discrepancies in upper bounds are not 
spurious. 

Being able to obtain minimum Hamming distance bounds 
as a function of the weight matrix is also interesting because 
the weight matrix is tightly connected to the concept of proto- 
graphs and LDPC codes derived from them 0, 0. Proto- 
graph-based code constructions start with a proto-graph that 
is described by a J x L incidence matrix whose entries are 
non-negative integers and where a "0" entry corresponds to 
no edge, a "1" entry corresponds to a single edge, a "2" entry 
corresponds to two parallel edges, etc. (Such an incidence 
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matrix is also known as proto-matrix.) Once such a proto- 
graph is specified, a proto-graph-based LDPC code is then 
defined to be the code whose Tanner graph iflOl is some r- 
fold graph cover ifTTl . Ifl2l of that proto-graph. 

It is clear that the construction of QC LDPC codes can then 
be seen as a special case of the proto-graph-based construction: 
first, the weight matrix corresponds to the proto-matrix, i.e., 
the incidence matrix of the proto-graph; secondly, the r-fold 
cover is obtained by restricting the edge permutations to be 
cyclic. 

A main reason for the attractiveness of QC LDPC codes is 
that they can be decoded using efficient decoding algorithms, 
namely belief-propagation-based decoding algorithms lfl3l or 
LP-based decoding algorithms OH, US), 02), El- Although 
the behavior of these decoders is mostly dominated by pseudo- 
codewords ED, 03, ED, ED, El, ED, El, E3 and 

the (channel-dependent) pseudo-weight of pseudo-codewords, 
the minimum Hamming distance still plays an important role 
because it characterizes undetectable errors and it provides 
an upper bound on the minimum pseudo-weight of a Tanner 
graph representing a code. 

The remainder of this paper is structured as follows. Sec- 
tion HI] introduces important concepts and the notation that will 
be used throughout the paper. Thereafter, Section [Til] presents 
the two main results of this paper. Both results are upper 
bounds on the minimum Hamming distance of a QC code: 
whereas in the case of Theorem[7]the upper bound is a function 
of the polynomial parity-check matrix of the QC code, in the 
case of Theorem [8] the upper bound is a function of the weight 
matrix of the QC code only. The following two sections are 
then devoted to the study of special cases of these results. 
Namely, Section [IV] focuses on so-called type-I QC LDPC 
codes (i.e., QC LDPC codes where the weight matrix entries 
are at most 1) and Section [V] focuses on so-called type-II 
and type-Ill QC LDPC codes (i.e., QC LDPC codes where 
the weight matrix entries are at most 2 and 3, respectively). 
We will show how we can obtain type-II and type-Ill codes 
from type-I codes having the same regularity and possibly 
better minimum Hamming distance properties. Section [VT] 
investigates the influence of cycles on minimum Hamming 
distance bounds. Finally, Section IVIII discusses a promising 
construction of type-I QC LDPC codes based on type-II or 
type-Ill QC LDPC codes. In fact, we suggest a sequence 
of constructions starting with a type-I code exhibiting good 
girth and minimum Hamming distance properties, or having 
good performance under message-passing iterative decoding. 
We construct a type-II or type-Ill code with the same regularity 
and higher Hamming distance upper bound, and from this we 
obtain a new type-I code with possibly larger minimum Ham- 
ming distance, up to twice the minimum Hamming distance 
of the type-II or type-Hi codes. Section IVIIII concludes the 
paper. The appendix contains many of the proofs and also one 
section (i.e, Appendix [J]) that lists some results with respect 
to graph covers. 

II. Definitions 

This section formally introduces the objects that were 
discussed in Section [J along with some other definitions that 



will be used throughout the paper. 

A. Sets, Rings, Fields, Vectors, and Matrices 

We use the following sets, rings, and fields: for any positive 
integer L, [L] denotes the set {0, 1, . . . , L— 1}; Z is the ring of 
integers; for any positive integer r, Z/rZ is the ring of integers 
modulo r; F2 is the Galois field of size 2; F2[x] is the ring of 
polynomials with coefficients in F2 and indeterminate x; and 
¥2[x}/ (x r — 1) is the ring of polynomials in ¥2[x] modulo x r — 
1, where r is a positive integer. We will often use the notational 
short-hand F 2 r '[x] for F 2 [x]/(x r -l). By F 2 \ F 2 r> [x] n , F™ xn , 
and Fj [x] mxn we will mean, respectively, a row vector over 
F2 of length n, a row vector over F 2 ^ [x] of length n, a matrix 
over F2 of size mxn, and a matrix over [x] of size m x n. 
In the following we will use the convention that indices of 
vector entries start at (and not at 1), with a similar convention 
for row and column indices of matrix entries. 

For any matrix M, we let .5 be the sub-matrix of M 
that contains only the rows of M whose index appears in the 
set 1Z and only the columns of M whose index appears in 
the set S. If 1Z equals the set of all row indices of M, we 
will omit in Mr s the set 1Z and we will simply write M5. 
Moreover, we will use the short-hand M^y; for M s \^j. 

As usual, the min operator gives back the minimum value 
of a list of valuesQ In the following, we will also use a more 
specialized minimum operator, namely the min* operator that 
gives back the minimum value of all nonzero entries in a list 
of values@ (Note that for lists containing only non-negative 
values, as will be the case in the remainder of this paper, the 
min* operator gives back the smallest positive entry of the 
list.) 

B. Weights 

The weight wt(c(x)) G Z of a polynomial c(x) G F 2 [x] 
equals the number of nonzero coefficients of c{x). Simi- 
larly, the weight wt(c(x)) G Z of a polynomial c(x) G 
F 2 r ^[x] equals the weight wt(c'(x)) of the (unique) minimal- 
degree polynomial c'(x) G F 2 [x] that fulfills c'(x) = 
c(x) (in F 2 r> [x]). 

Let c(x) = (co(x),cx(x),...,c L _x(x)) G F 2 r> [x] L 
be a length-!/ polynomial vector. Then the weight vector 
wt(c(x)) G Z L of c(x) is a length-i vector with the i-th 
entry equal to wt(ci(x)). Similarly, let H(x) = \hj : i(x)\. . G 

F 2 r ^[x] JxL be a size-JxL polynomial matrix. Then the 
weight matrix wt(H(x)) G Z ,7xL of H(x) is a JxL-matrix 
with the entry in row j and column i equal to wt[/i Jj j(x)] . 

The Hamming weight wh(c) of a vector c is the number 
of nonzero entries of c. In the case of a polynomial vector 
c(x) = (co(x), ci(x), . . . , Cl-i(x)) G F 2 [x] l , the Hamming 
weight wh(c(x)) is defined to be the sum of the weights of 
its polynomial entries, i.e., wh(c(x)) = J2f=a ( w ^( c ( x ))) ■ = 
T,f=o wt(cj(x)). 

'If the list is empty then min gives back +00. 

2 If zero is the only value appearing in the list then min* gives back +00. 
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C. Codes 

All codes in this paper will be binary linear codes. As usual, 
a code C of length n can be specified through a (scalar) parity- 



check matrix H 6 



jpm X n 



i.e. 



C = {c e F5 



H 



c = 



T }, where T denotes transposition. This code has rate at least 
1 — — and its minimum Hamming distance (which equals the 
minimum Hamming weight since the code is linear) will be 
denoted by e? min (C). 

Let J, L, and r be positive integers. Let C be a code of 
length Lr that possesses a parity-check matrix H of the form 



H 



H 
H 



o.o 

1.0 



H .i 
Hii 



Ho,L-l 

Hi i_i 



where the sub-matrices H 



H 



J-l.L-l 



e F 



JrxLr 



j i E F" 1, are circulant. Such a code 



is called quasi-cyclic (QC) because applying circular shifts to 
length-r subblocks of a codeword gives a codeword again. 
Because H, $ is circulant, it can be written as the sum H 



si 



-0 ^M,s,0 



I s , where ft.j,i, s ,o is the entry of Hj 



J; 1 ~ 

in row 



s and column 0, and where I s is the s times cyclically left- 
shifted identity matrix of size r x r. 

With a parity-check matrix H e F 2 7rxir of a QC code 
we associate the polynomial parity-check matrix H(.t) G 



F 2 r> [x] JxL 



H(aO 



ho,o( x ) 
hi,o(x) 



hiUx) 



hi,L-i(x) 



hj- 1>0 (x) hj-i tl (x) 
where h jA (x) = YSsZl hj ti , sfi x s . 

c = (C0,0j • • • ) C0,r-lj • • • i c i-l,0j ■ 

associate the polynomial vector 



hj. 



1,L-1 



(X\ 



Moreover, to any vector 

..,C£_i, r _i) e ¥% r we 



c(z) = (cd(s),ci(x),... )C i_i(x)) eff[x]" , 

where Ci(x) = 5Zs=o °i,sX s - K can easily be checked that the 
condition 

H • c T = T (in F 2 ) 
is equivalent to the condition 

H(x) ■ c(x) T = T (in Fi, r> [x]) , 



giving us an alternate way to check if a (polynomial) vector 
is a codeword. 

The following classification was first introduced in (7). 

Definition 1. Let M be some positive integer. We say that a 
polynomial parity-check matrix H(x) of a QC LDPC code is 
of type M if all the entries of the associated weight matrix 
wt(H(a:)J are at most M. □ 

Equivalently, H(x) is of type M if for each polynomial 
entry in H(x) the number or nonzero coefficients is at most 
M. In particular, the polynomial parity-check matrix H(x) 
is of type 1 (in the following often called type I) if H(x) 
contains only the zero polynomial and monomials. Moreover, 
the polynomial parity-check matrix H(x) is of type 2 (in the 



following often called type II) if H(x) contains only the zero 
polynomial, monomials, and binomials. If H(x) contains only 
monomials then it will be called a monomial parity-check 
matrix. (Obviously, a monomial parity-check matrix is a type-I 
polynomial parity-check matrix.) 

D. Graphs 

With a parity-check matrix H we associate a Tanner 
graph IfTUl in the usual way, namely, for every code bit we 
draw a variable node, for every parity-check we draw a check 
node, and we connect a variable node and a check node by an 
edge if and only if the corresponding entry in H is nonzero. 
Similarly, the Tanner graph associated to a polynomial parity- 
check matrix H(x) is simply the Tanner graph associated to 
the corresponding (scalar) parity-check matrix H. 

As usual, the degree of a vertex is the number of edges 
incident to it and an LDPC code is called (d\, e^) -regular if 
all variable nodes have degree di and all check nodes have 
degree c^. Otherwise we will say that the code is irregular. 
Moreover, a simple cycle of a graph will be a backtrackless, 
tailless, closed walk in the graph, and the length of such a 
cycle is defined to be equal to the number of visited vertices 
(or, equivalently, the number of visited edges). The girth of 
a graph is then the length of the shortest simple cycle of the 
graph. 

The above-mentioned concepts are made more concrete with 
the help of the following example. 

Example 2. Let C be a length- 12 QC LDPC code that is 
described by the parity-check matrix 
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Clearly, J = 3, L = 4, and r = 3 for this code and so H can 
also be written like 



H 



Io + Ii Io 
I 2 Io 
Ii 



I 2 
Ii I 2 
Io + h Ii 



where I s , s = 0, 1, . . . , r— 1, are s-times cyclically left-shifted 
r x r identity matrices. The corresponding polynomial parity- 
check matrix is 



H(x) = 
and the weight matrix is 



x + x 

^2 







x" + x 2 



wt (H(a?)) = 



2 10 1 
1111 
12 1 
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Fig. 1. Left: Tanner graph of a length-12 QC LDPC code. It is a triple cover 
of the proto-graph shown on the right. Right: Proto-graph of the Tanner graph 
shown on the left. 

The Tanner graph associated to H or H(x) is shown in 
Figure Q] (left). We observe that all variable nodes have degree 
3 and all check nodes have degree 4, therefore C is a (3, 4)- 
regular LDPC code. (Equivalently, all columns of H have 
weight 3 and all rows of H have weight 4.) □ 

The proto-graph associated to a polynomial parity-check 
matrix H(x) is a graphical representation of the weight matrix 
wt(H(x)) in the following way. It is a graph where for each 
column of H(a;) we draw a variable node, for each row of 
H(x) we draw a check node, and the number of edges between 
a variable node and a check node equals the corresponding 
entry in wt(H(x)). 

Example 3. Continuing Example 12 the proto-graph of H(.t) 
is shown in Figure [T] (right). Clearly, the weight matrix 
wt (H(a;)) is the incidence matrix of this latter graph. We 
observe that all variable nodes have degree 3 and all check 
nodes have degree 4. (Equivalently, all column sums (in Z) 
of wt (H(a:)) equal 3 and all row sums (in Z) of wt (H(x)) 
equal 4.) □ 

An important concept for this paper are so-called graph 
covers, see the next definition. 

Definition 4 (See, e.g., IfTTI . |[T2l ). Let G be a graph with 
vertex set V(G) and edge set £ (G), and let d(v) denote the 
set of adjacent vertices of a vertex v G V(G). An unramified, 
finite cover, or, simply, a cover of a (base) graph G is a graph 
G along with a surjective map (f> : G — ► G, which is a graph 
homomorphism, i.e., which takes adjacent vertices of G to 
adjacent vertices of G such that, for each vertex v G V(G) 
and each v G the neighborhood d(v) ofv is mapped 

bijectively to d(v). For a positive integer r, an r-cover of G is 
an unramified finite cover <p '■ G — > G such that, for each vertex 
v G V(G) of G, <j> (v) contains exactly r vertices of G. An 
r-cover of G is sometimes also called an r-sheeted covering 
of G or a cover of G of degree 

3 It is important not to confuse the degree of a covering and the degree of 
a vertex. 



Example 5. Continuing Examples [2] and [3] we note that 
the graph in Figure [T] (left) is a 3-cover of the graph in 
Figure [TJ (right). Therefore, the code C is a proto-graph-based 
code. It can easily be checked visually that all the edge 
permutations that were used to define this 3-cover are cyclic 
permutations, therefore confirming that the code C is indeed 
quasi-cyclic. □ 

We conclude this subsection by emphasizing that graph 
covers have been used in two different ways in the context 
of LDPC codes: on the one hand, they have been used for 
constructing LDPC codes (like in this paper), on the other 
hand they have been used to analyze message-passing iterative 
decoders (like in ET1. El). 

E. Determinants and Permanents 

The determinant of an m x m-matrix B ~ over 
some ring is defined to be 

det(B)=£(-lf^) J] b jMj) , 

where the summation is over all ml permutations of the 
set [m], and where (— l) sl s n (< 7 ) equals +1 if a is an even 
permutation and equals — 1 if a is an odd permutation. 

The permanent of an m x ?n-matrix B = \bj,i]j.i over some 
ring is defined to be 

perm(B) = ^ II b i,*(j) > 

a je[m] 

where the summation is over all ml permutations of the set 

[m\. 

Clearly, for any matrix B with elements from a ring or field 
of characteristic 2 it holds that dct(B) = pcrm(B). 

III. Minimum Hamming Distance Upper Bounds 

This section contains the two main results of this paper, 
namely Theorems [7j and [8] More precisely, given some quasi- 
cyclic LDPC code with polynomial parity-check matrix H(.t) 
and minimum Hamming distance d m i n (C), Theorem[7jpresents 
an upper bound on <i m i n (C) as a function of the entries of 
H(x) and Theorem [8] presents an upper bound on d m i n (C) 
as a function of the entries of wt(H(x)). The upper bound 
in Theorem [8] is in general weaker than the upper bound 
in Theorem [7] however, it is interesting to see that the 
weight matrix alone can already give nontrivial bounds on 
the achievable minimum Hamming distance. 

In later sections, namely in Sections |IV] and [V] we will 
discuss the implications of these two theorems on codes with 
type-I, type-II, and type-Ill polynomial parity-check matrices. 
Moreover, in Section [VI] we will show how the upper bounds 
in Theorems[7jand[8]can be strengthened by taking some graph 
structure information (like cycles) into account. 

We start with a simple technique presented by MacKay 
and Davey |4) to construct codewords of codes described 
by polynomial parity-check matrices. (Actually, the paper (4) 
deals with codes that are described by scalar parity-check 
matrices composed of commuting permutation submatrices, of 
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which parity-check matrices composed of cyclically shifted 
identity matrices are a special case. However, and as we 
show in this paper, their techniques can be suitably extended 
to codes that are described by scalar parity-check matrices 
composed of any circulant matrices, and therefore to codes 
that are described by polynomial parity-check matrices.) 

Lemma 6. Let C be a QC code with a polynomial parity- 
check matrix H(x) G F% [x] JxL whose entries are de- 
noted by hj.i(x), i.e., H(x) = [hj,i{%)]j,i- Let <S be 
an arbitrary size-(J+l) subset of [L] and let c(x) = 
(co(x), Ci(x), . . . , cl-i(x)) G [x]^ be a length-L vector 
defined byj 

(^(x) = |P crm ( H sv( x )) i f i e S 
1 otherwise 

Then c(x) is a codeword in C. 

Proof: See Appendix lAl ■ 

With the help of the codeword construction technique in 
Lemma|6]we can easily obtain the bound in Theorem[7] simply 
construct the list of all codewords corresponding to all size- 
(J+l) subsets S of [L], and use the fact that the minimum 
Hamming distance of C is upper bounded by the minimum 
Hamming weight of all non-zero codewords in this list. 

Theorem 7. Let C be a QC code with polynomial parity-check 
matrix H(x) <E F^ r> [x] JxL . Then the minimum Hamming 
distance of C is upper bounded as follows 

dmin(C) < min* Y^wt (perm(H 5 \ l (.T)) ) . (1) 

\s\=.i+i ies 

Proof: Let S be a size- (J +1) subset of [L] and let 
c(x) be the corresponding codeword constructed according 
to Lemma [6] The result in the theorem statement follows by 
noting that c(x) has Hamming weight 

wn(c(x)) = ^2 wt(a(x)) = 2jwt (ci(x)) 
ie[L] ies 

= 5Z wt (p erm ( H 5\i( a; ))) ■ 



Let us emphasize that it is important to have the min* 
operator in (fl]i, and not just the min operator. The reason 
is that the upper bound is based on constructing codewords of 
the code C and evaluating their Hamming weight. For some 
polynomial parity-check matrices some of these constructed 
codewords may equal the all-zero codeword and therefore have 
Hamming weight zero: clearly, such constructed codewords are 
irrelevant for upper bounding the minimum Hamming distance 
and therefore must be discarded. This is done with the help 
of the min* operator. 

The next theorem, Theorem [8] gives an upper bound on 
bound (03 which is easier to compute and depends only on 

4 Because the ring [x] has characteristic 2, we could equally well define 
Ci{x) = det (H^^fi)) if % 6 S. 



the weight matrix associated to H(x). In particular, this bound 
does not depend on r, the size of the circulant matrices in 
the scalar parity-check matrix H corresponding to H(x). The 
bound says that the minimum Hamming distance is upper 
bounded by the minimum nonzero sum of the permanents of 
all J x J submatrices of a chosen J x (J+l) submatrix of 
the weight matrix A, the minimum being taken over all such 
possible J x (J+l) submatrices of A. 

Theorem 8. Let C be a QC code with polynomial parity-check 
matrix H(x) G F^[x] JxL and let A = wt(H(a;)). Then the 
minimum Hamming distance of C is upper bounded as follows 

drain (C) < nun* ^ P erm ( A 5\i) ■ (2) 

•si=j+i ies 

Proof: See Appendix 151 ■ 

Again, as in Theorem [7] it is important to have the min* 
operator in Theorem[8]and not just the min operator. This time 
the reasoning is a bit more involved, though, and we refer the 
reader to the proof of Theorem [8] for details Jf| 

Note that the upper bound in (Q]) depends on r (because the 
computations are done modulo X r — 1), whereas the bound 
in (ffj does not depend on r. In fact, the latter bound is a bound 
on the free Hamming distance of the convolutional code that 
could be associated to this family of QC codes, which is an 
upper bound on the minimum Hamming distances of the QC 
codes in this family, see J27), EH). 

Based on Theorems [7] and [8] the following recipe can 
be formulated for the construction of QC codes with good 
minimum Hamming distance. 

• Search for a suitable weight matrix with the help of 
Theorem [8] 

• Among all polynomial parity-check matrices with this 
weight matrix, find a suitable polynomial parity-check 
matrix with the help of Theorem UJ 

• Verify explicitly if the minimum Hamming distance of 
the code of the found polynomial parity-check matrix 
really equals (or comes close to) the minimum Hamming 
distance promised by the upper bound in Theorem UJ 

This recipe is especially helpful in the case where one is 
searching among type-M polynomial parity-check matrices 
with small M, say M £ {1,2,3}. In such cases it is to 
be expected that there is not much difference in the upper 
bounds (Q]i and (0. For type-M polynomial parity-check 
matrices with larger M, however, we do not expect that the 
upper bounds (Q~|) and (f2]) are close. The reason is that when 
computing permfH^x i(x)) in (UJ there will be many terms 
that cancel each other. Anyway, when constructing QC LDPC 
codes, type- M polynomial parity-check matrices with large 
M are somewhat undesirable because of the relatively small 
girth of the corresponding Tanner graph. In particular, it is 
well known that for M 3 the Tanner graph of a polynomial 
parity-check matrix whose weight matrix contains at least one 
entry of weight M has girth at most 6 (see also Theorem [171). 

5 We are grateful to O. Y. Takeshita for pointing out to us that in earlier 
(and also less general) versions of Theorem [7] and Theorem [8] (cf. 0) the 
min operator has to be replaced by the min* operator, see also 1261 . 
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IV. Type-I QC Codes 

In this section we specialize the results of the previous 
section to the case of type-I parity-check matrices. 

Corollary 9. Let C be a type-I QC code with polynomial 



WuJxi 



and let A 



parity-check matrix H(x) £ ¥2 
wt(H(ai)). Then the minimum Hamming distance ofC is upper 
bounded as follows 



d min (C) < (J+l)! 
Proof: See Appendix ICl 



(3) 



Let us evaluate the minimum Hamming distance upper 
bounds that we have obtained so far for some type-I poly- 
nomial parity-check matrices. (Actually, the following poly- 
nomial parity-check matrices happen to be monomial parity- 
check matrices, i.e., polynomial parity-check matrix where all 
entries of the corresponding weight matrix equal 1.) 

Example 10. Let r ^ 15. Consider the (2, 4)-regular length- 
4r QC-LDPC code C given by the polynomial parity-check 
matrix H(x) 



H(x) 



x 

„5 



and the (2, 4)-regular length-4r QC-LDPC code C given by 
the the polynomial parity-check matrix 



H'( 



x 
>,<; 



According to (Q}, the minimum Hamming distance of C is 
upper bounded by 

'wt(a; 7 +a; 7 ) + wt(a; 5 +x 10 ) + wt(x 4 +x 9 ), 
■ * J wt(ir 7 +x 7 ) + wt(x 8 +x 13 ) + wt(x 9 +x 14 ), 
rf min ^ mm < wt(a ,4 +s 9) + ^(^8+3^3) + wt(.T 11 +a; 11 ), 

Wt{x 5 +X 10 ) + Wt(.T 9 +X 14 ) + Wt(x 11 +.T 11 ) j 

= min*{4,4,4,4} = 4 , 

and the minimum Hamming distance of C is upper bounded 
by 

{wt(x 6 +.T 8 ) +wt(x 5 +x 9 ) + wt(a; 4 +a; 10 ), 
wt(x 6 +x s ) + wt(x 10 +x 14 ) +wt(x 11 +x 13 ), 
wt(x 4 +x w ) + wt(x w +x 14 ) + wt(x 13 +z u ) 
wt(x 5 +x 9 ) + wt(a- n +a; 13 ) + wt(x 13 +x n ) 

= min{6, 6, 6, 6} = 6 . 
However, in both cases the bound in (f2]i gives 

^min 

{(1+1) + (1+1) + (1+1)} =6 

since both polynomial parity-check matrices have the same 
weight matrix. Similarly, in both cases the bound in (O gives 

dmin < (2 + 1)! = 6 

since both polynomial parity-check matrices have J = 2. 

In conclusion, we see that a 2x4 monomial parity-check 
matrix can yield a code with minimum Hamming distance at 
most 6. However, when the entries of the polynomial parity- 
check matrix are not chosen suitably, as is the case for H(x), 



then the minimum Hamming distance upper bound in (HJ is 
strictly smaller than the minimum Hamming distance upper 
bound in (fj}. 

For completeness, we computed the minimum Hamming 
distance of the two codes, and obtained 2 for the first code, as 
(0, 0, x 4 , 1) is a codeword, and 4 for the second code for most 
values of r (the free Hamming distance of the convolutional 
code with parity-check polynomial matrix H'(D) is 4). □ 

Let us discuss another example. 

Example 11. Let r ^ 21 and let the (3, 4) -regular QC- 
LDPC code C be given by the polynomial parity-check matrix 

H(a:) G : 



13x4 



X 


X 2 


X 4 


x s 


X 5 


x w 


x 20 


X 9 


X 25 


x 19 


x 7 


x 14 



H(x) 



(This code was obtained by shortening the last r positions 
of the (3, 5)-regular type-I QC-LDPC code of length 5r 
presented in [;29||@) Evaluating the bounds in (ffjl and §3^ for 
this polynomial parity-check matrix, we see that the minimum 
Hamming distance is upper bounded by 24, and for suitable 
choices of r this upper bound is indeed achieved. We computed 
the minimum Hamming distance of the code for different 
values of r and obtained that r = 31 is the smallest such 
choice. The minimum Hamming distance for some values are 
listed in the following table. 



r 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


drain ) 


10 


8 


20 


16 


18 


18 


14 


16 


18 


8 


24 



□ 

As we have seen from the above examples, the minimum 
Hamming distance upper bound (Q]i can be strictly smaller than 
the upper bound However, the latter minimum Hamming 
distance upper bound is computed more easily, and it provides 
an upper bound on the Hamming distance of all QC codes 
having the same weight matrix and therefore also the same 
basic proto-graph. 

Applying Corollary [9] to codes with monomial parity-check 
matrices shows that for such codes the minimum Hamming 
distance is upper bounded by (J+l)!. We note that this result 
was previously presented by MacKay and Davey 0] and 
discussed by Fossorier J5). (As already mentioned before, 
MacKay and Davey [4 | focused on codes that are described by 
scalar parity-check matrices composed of commuting permuta- 
tion submatrices, of which parity-check matrices composed of 
cyclically shifted identity matrices are a special case. However, 
and as we show in this paper, their techniques can be suitably 
extended to codes that are described by scalar parity-check 
matrices composed of any circulant matrices, and therefore to 
codes that are described by polynomial parity-check matrices.) 

V. Type-II and Type-III QC Codes 

After having discussed minimum Hamming distance upper 
bounds for type-I QC codes in the previous section, we now 

6 Note that in 1291 , r = 31, and the code parameters are [155, 64, 20]. Also 
note that by shortening a code, the girth of the associated Tanner graph cannot 
decrease. 
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present similar results for type-II and type-Ill QC codes. In 
particular, we classify all possible weight matrices of (3,4)- 
regular QC codes with a 3 x 4 polynomial parity-check matrix. 

We start our investigations with the following motivating 
example. 

Example 12. In Example QT| we saw that the minimum 
Hamming distance of type-I (3, 4)-regular QC codes with a 
3x4 polynomial parity-check matrix cannot surpass 24. In this 
example we show that type-II (3, 4)-regular QC codes with a 
3x4 polynomial parity-check matrices can have minimum 
Hamming distance strictly larger than 24. Namely, consider 
the code C with parity-check matrix 



U'(x) 



(This polynomial parity-check matrix was obtained from the 
parity-check matrix H(x) in Example [TT| by pairing some 
monomials into binomials and replacing with the positions 
left, careful to preserve the (3, 4)-regularity.) The correspond- 
ing weight matrix is 



upper bound implied by (01 ): 



+ x 2 





x i 


X 8 




X 5 


X 9 


x w + x 20 





■ (4) 





x 25 + x 19 





X 7 + x 14 





A' = 



2 11 
112 
2 2 



and, according to (0, yields the following minimum Hamming 
distance upper bound 



d m m (C) < mm* {10 + 6 + 10 + 6} = 32 . 



For small r, the corresponding code does not attain this bound, 
however, for r = 46 one can verify that the resulting code 
attains the optimal minimum Hamming distance d m j n = 32. 
□ 

After this introductory example, let us have a more system- 
atic view of the possible weight matrices of (3, 4)-regular QC 
codes and the minimum Hamming distance upper bounds that 
they yield. 

Corollary 13. Let C be a (3, 4)-regular type-II QC code with 
polynomial parity-check matrix H(x) G F^fx] 3 * 4 and let 
A = wt(H(x)) G Z 3x4 . Then all possible (3, A)-regular type- 
II weight matrices A. (up to permutations of rows and columns) 
are given by the following 5 types of matrices ( shown here 
along with the corresponding minimum Hamming distance 



2 2 
1111 




2 2 



2 2 





1 1 




2 1 



2 

1 1' 



1 1 
1 1 



1111 

1111 
1111 
1111 



with d-n 



8 + 8 =32 



with d min < 10 + 10 + 6 + 6 = 32 , 



with d mi „ ^7+7+7+9 =30, 



with d m i n 5:6 + 6 + 8 



= 28, 



with d min ^6 + 6 + 6 + 6 =24. 



As can be seen from this list, the largest minimum Hamming 
distance upper bound is d m in ^ 32 and it can be obtained 
if the weight matrix A equals ( modulo permutations of rows 
and columns) the first or the second matrix in the list. 

Proof: Omitted. ■ 

Example 14. We see that the type-II (3, 4)-regular QC code 
with a 3 x 4 polynomial parity-check matrix and with r = 46 
presented in Example[l2lnot only achieves the minimum Ham- 
ming distance upper bound promised by (f2]i but, according to 
Corollary [13] it achieves the best possible minimum Hamming 
distance upper bound for any type-II (3, 4)-regular QC code 
with a 3 x 4 polynomial parity-check matrix. We note that 
this particular code has parameters [184, 47, 32], girth 8, and 
diameter 8, i.e., the same girth and diameter as the Tanner 
graph of the [124, 33, 24] code in Example fTTb . which is a 
shortened version of the [155,64,20] code in 1291 . 

Figure [2] shows the decoding performance (word error 
rate) of the [184,47,32] QC LDPC code when used for 
transmission over a binary-input additive white Gaussian noise 
channel. Decoding is done using the standard sum-product 
algorithm fPUl which is terminated if the syndrome of the 
codeword estimate is zero or if a maximal number of 64 
(respectively 256) iterations is reached. It is compared with a 
randomly generated (3, 4)-regular [184, 46] LDPC code where 
four-cycles in the Tanner graph have been eliminated. (When 
comparing these two codes one has to keep in mind that 
because the randomly generated code has slightly lower rate 
and because the horizontal axis shows E-^/Nq, the randomly 
generated code has a slight "disadvantage" of 0.093 dB.) Note 
though that the decoding complexity per iteration is the same 
for both codes. 

Let us mention on the side that we tried to estimate the 
minimum (AWGN channel) pseudo-weight 11211 . Il22l of this 
code. From searching in the fundamental cone we get an 
upper bound of 27.6 on the minimum pseudo-weight. The 
pseudo-weight spectrum gap l24l is therefore estimated to be 
32 — 27.6 = 4.4, which is on the same order as the pseudo- 
weight spectrum gap for the (3, 5)-regular [155, 64, 20] code 
by Tanner |29l , which is estimated to be 20 — 16.4 = 3.6. We 
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e Random code (after max. 64 iterations) 
Random code (after max. 256 iterations) 
A QC-LDPC code (after max. 64 iterations) 
-A- QC-LDPC code (after max. 256 iterations) 



10' 



E b /N [dB] 



Fig. 2. Decoding performance of the [184, 47, 32] QC-LDPC code vs. a randomly generated (four-cycle free) [184, 46] LDPC code under sum-product 
algorithm decoding when transmitting over a binary-input AWGN channel. (For more details, see Example 1 14h 



also note that for the above-mentioned randomly generated 
[184, 46] code we obtained an upper bound on the minimum 
pseudo-weight of 21.0. □ 

If we want to take into consideration all cases of (3,4)- 
regular QC codes, we also have to investigate the class of 
type-Ill QC codes, as is done in the next corollary. 



matrix A equals (modulo row and column permutations) the 
last matrix in the above list. 

Proof: Omitted. ■ 

Example 16. We can modify the matrix H in Example [Til to 
obtain one of the configurations in Corollary Q3] For example, 



the following matrix H'(x) G F 



(r). 



13x4 



Corollary 15. Let C be a (3, 4) -regular type-Ill QC code last configuration listed in Corollary! 

with polynomial parity-check matrix H(x) G F, 



corresponds to the 



13x4 



and 



let A = wt(H(ai)) G Z 3x4 . Then all possible (3,4)-regular 
type-Ill weight matrices A (up to permutations of rows and 
columns) are given by the 5 types of matrices already listed in 
Corollary\T5\ together with the 3 types of matrices (shown here 
along with the corresponding minimum Hamming distance 
upper bound implied by ©J: 



H'(x) 4 



-x 







,10. 





x 

X 5 



3 

2 

1 

'3 

2 





1 

1 1 

2 1 

0' 

1 1 

2 2 



3 1 
3 1 
3 1 



with d min <5 + 9 + 9 + 15 =38 



with dmin < 4 + 12 + 12 + 12 = 40 



with d min ^9 + 9 + 9 + 27 =54 



For r = 31 we obtain a [124, 31, 28] code, whose rate is 0.25. 
(In comparison, the monomial [124, 33, 24] QC LDPC code 
in Example [UJ has rate 0.266 and the binomial [184,47,32] 
QC LDPC code in Example [TJ has rate 0.2554.) For r = 46, 
we obtain a trinomial code with parameters [184,46,34]. For 
larger r the minimum Hamming distance could increase up to 
54. □ 

Note that the Tanner graph of a polynomial parity-check 
matrix which has at least one trinomial entry cannot have girth 
larger than 6 (see, e.g., 0). We state this observation as part 
of a more general analysis of the effect of the weight matrix 
on the girth. 



As can easily be seen, the largest minimum Hamming distance 
upper bound is d m [ n ^ 54 and it can be obtained if the weight 



Theorem 17. Let C be a QC code described by a polynomial 
parity-check matrix H(x) G [x] Jxi . Let gir be the girth 
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of the Tanner graph corresponding to H(x) and let A be the to the girth of their Tanner graph: firstly, there is a lot of 



weight matrix corresponding to H(x). 



a) If A. has submatrix 

b) If A. has submatrix 



1 1 1 
1 1 1 



1 1 
1 2 



then gir ^12. 
then gir ^10. 



c) If A has submatrix [2 2] then gir ^ 8. 

A) If A. has submatrix [3] then gir 6. 
(By "A having submatrix B" we mean that A contains a 
submatrix that is equivalent to B, modulo row permutations, 
column permutations, and transposition.) 

Proof: See Appendix ID1 ■ 

The Corollaries [131 and [151 focused on the case of (3,4)- 
regular QC codes with a 3 x 4 polynomial parity-check matrix. 
It is clear that similar results can be formulated for any 
(</', Z7)-regular QC code with a JxL polynomial parity-check 
matrix. However, we will not elaborate this any further, except 
for mentioning the following corollary about (3, 5)-regular QC 
codes with a 3 x 5 polynomial parity-check matrix. 

Corollary 18. An optimal (3, b)-regular type-II weight matrix 
of size 3x5 must (up to row and column permutations) look 
like 

"2 2 1 0" 
A= 2 2 1 . 
110 12 

This weight matrix yields the minimum Hamming distance 
upper bound 

d m in «S min{30, 30, 30, 32, 28} = 28 . 

Proof: Omitted. ■ 



VI. The Effect of Small Cycles 
on the Minimum Hamming Distance 

If we know that the Tanner graph corresponding to some 
polynomial parity-check matrix contains some short cycles 
then we can strengthen the upper bounds of Theorem [7] and [8] 
In particular, Theorems [21] [24] and [25] will study the influ- 
ence of 4-cycles, 6-cycles, and 2i?-cycles, respectively, upon 
minimum Hamming distance upper bounds. These theorems 
will be based on results presented in Lemmas [19] and [23j 
that characterize cycles in Tanner graphs in terms of some 
entries of the corresponding polynomial parity-check matrix, 
especially in terms of permanents of submatrices. In order to 
state such conditions, we will use results from ||29l . 0. (For 
other cycle-characterizing techniques and results, see also ll30l 
and ED.) 

As we will see, the smaller the girth of the Tanner graph, 
the smaller the minimum Hamming distance upper bound 
will be. This observation points in the same direction as 
other results that relate the decoding performance of LDPC 
codes, especially under message-passing iterative decoding, 



empirical evidence that smaller girth usually hurts the iterative 
decoding performance; secondly, there are results concerning 
the structure of the fundamental polytope that show that the 
fundamental polytope of Tanner graphs with smaller girth is 
"weaker" 



A. Type-I QC Codes with 4-Cycle 

Lemma 19. The Tanner graph of a type-I QC code with 
polynomial parity-check matrix H(x) € F^ [x] JxL has a 4- 
cycle if and only ifH(x) has a 2x2-submatrix B(x)/or which 



wt ^pcrm (B(a;))^ < perm ( wt (B(a;)) 



holds. 



Proof: See Appendix [E] 



Corollary 20. The Tanner graph of a type-I QC code with 



Wr 



[JxL 



has a 4- 



polynomial parity-check matrix H(x) G F 
cycle if and only ifH(x) has a 2x2-submatrix B(a;) which is 
monomial and for which perm(B(a;)) ~ (in ¥^[x]) holds. 

Proof: Follows from Lemma [19] and its proof. ■ 

With this, we are ready to investigate the minimum Ham- 
ming distance upper bound for Tanner graphs with 4-cycles. 

Theorem 21. Let C be a type-I QC code with polynomial 



parity-check matrix H(x) G F. 



<r>. 



1 JxL 



If the associated 



Tanner graph has a A-cycle then 

draUQ < (•/+!)! - 2(J-1)! . 
Proof: See Appendix [F] 



(5) 



Example 22. Let us consider again the (2, 4) -regular length- 
4r QC code C from Example [10] which is given by the 
polynomial parity-check matrix 

H(x) 



It has two 4-cycles since 



perm 



= and perm 



= 



(in ¥ { 2 r> [x}). Therefore the bound © gives 

dmin(C) ^ (J+l)! - 2(J-1)! = 3! - 2 ■ 1! 



(6) 



We note that for this H(x) this upper bound equals the upper 
bound © (cf. Example \TU$. □ 



B. Type-I QC Codes with 6-Cycles 

Lemma 23. The Tanner graph of a type-I QC-LDPC code 



1 JxL 



has 



with polynomial parity-check matrix H(ir) G F 
a 6-cycle (or possibly a 4-cycle) if and only if H(x) has a 
si'ze-3x3 submatrix B(x) for which 



wt 



perm (B(x)) ) < perm I wt (B(x)) 
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holds, i.e., if and only ifH(x) has a 3x3-submatrix B(x)/or 
which some terms of its permanent expansion add to zero. 

Proof: See Appendix iGl ■ 

With this, we are ready to investigate the minimum Ham- 
ming distance upper bound for Tanner graphs with 6-cycles. 

Theorem 24. Let C be a type-I QC code with polynomial 
parity-check matrix H(x) £ [x] JxL . If the associated 
Tanner graph has a 6-cycle then 

dmin < (J+l)! - 2(J-2)! . 
Proof: See Appendix IH1 ■ 



C. Type-I QC Codes with 2 J' '-Cycles 

The previous two subsections have shown that the minimum 
Hamming distance of a type-I code whose Tanner graph has 
girth 4 or 6 can never attain the maximal value (J+l)! of 
Corollary [9] These results are special cases of a more general 
result that we will discuss next. Note however that, compared 
to the girth-4 and the girth-6 case, this more general statement 
is uni-directional. 

Theorem 25. Let C be a type-I QC code with polynomial 
parity check matrix H(x) = [hj,i(x)]i t j € ¥^[x] JxL . Let R, 
2 ^ R ^ min( J, L), be some integer and suppose there is 
a set 1Z C [J] of size R, a set S C [L] of size R, and two 
distinct bijective mappings a and r from 1Z to S such that 
a {j) 7^ T (j) f or Q tt j £ ^ and such that 

n h j,°{j)( x ) = n h J,r(j)( x ) ■ ( ? ) 

Then 

cUm < (J+l)! - 2{J-R+l)\ . 

If, in addition, the (bijective) mapping a^ 1 or from 1Z to 1Z is 
a cycle permutation of order R and if the products on the left- 
hand and right-hand side of the equation in (O are nonzero, 
then the associated Tanner graph will have a cycle of length 
2R. 

Proof: See Appendix U ■ 

For 4- and 6-cycles the converse of the second part of the 
above corollary is true, i.e., 4- and 6-cycles are visible in, 
respectively, 2x2- and 3 x 3-submatrices (cf. Theorems |2TI 
andl24l. However, for longer cycles the converse of the second 
part of the above corollary is not always true: 8-cycles can 
happen in 4 x 4-submatrices, but also in 2 x 4-submatrices 
or in 3 x 4-submatrices. A similar statement holds for longer 
cycles. 

D. Type-II QC Codes 

Similar statements as in the preceding subsection can also 
be made about type-II QC LDPC codes. We will not say 
much about this topic except for stating a lemma that helps in 
detecting if a polynomial parity-check matrix is 4-cycle free. 



Lemma 26. A type-II QC code is 4-cycle free if and only 
if the polynomial parity-check matrix associated to it has the 
following properties. 

1) If r is even, then for any Ixl-submatrix like 

[x a +x b ] , 

it holds that 

gcd(x a + x b 1 l + x r ) ^ 1 + x r ' 2 (in ¥ 2 [x}) . 

2) For any 1 x 2-submatrix like 

[x a + x b x c + x d ] , 

the product (x a + x b ) ■ (x c + x d ) (in F^ [x]) has weight 
4, i.e., the maximally possible weight. 

3) For any 2x\-submatrix like 



x v 



the product (x a + x b ) ■ (x c + x d ) (in F^ [x]) has weight 
4, i.e., the maximally possible weight. 
4) For any 2 x 2-submatrix like 



x d + x e 



X + X 

T.f 



(or row and column permutations thereof), the perma- 
nents (in F^ [x]) of the following two 2 x 2-submatrices 

„b~\ 

and 



r.f 



have weight 3, the maximally possible weight. 
5) For any 2 x 2-submatrix with weight matrix 



2 


2 




'2 1" 


, and 


'l 1" 




1 




1 




1 1 


1 1 



(or row and column permutations thereof), the perma- 
nent (in F2 [x]) of this 2 x 2-submatrix has weight 4, 3, 
and 2, respectively, i.e., the maximally possible weight. 

Proof: It is well known that a 4-cycle appears in a Tanner 
graph if and only if the corresponding (scalar) parity-check 
matrix contains the 2x2 (scalar) submatrix 



The lemma is then proved by studying all possible cases in 
which a polynomial parity-check matrix can lead to a (scalar) 
parity-check that contains this this 2x2 (scalar) submatrix. 
The details are omitted. ■ 

Example 27. Consider the type-II polynomial parity-check 
matrix H'(x) in ©. There, the 4-cycles could only appear 
from the two submatrices 



x 



x 



x 



„10 



,.20 



and [a 



25 



,.10 



x x + x 

Therefore, the conditions for the non-existence of a 4-cycle 
are 

$ {5 - 1, 5 - 2} + {4 - 10, 4 - 20} (in Z/rZ) , 

i {25 - 19, 19 - 25} + {7 - 14, 14 - 7} (in Z/rZ) . 
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(Here the sum of two sets denotes the set of all possible sums 
involving one summand from the first set and one summand 
from the second set.) It is clear that similar analyses could be 
made for the non-existence of longer cycles. □ 

VII. Type-I QC Codes based on 
Double Covers of Type-II QC Codes 

So far, we have mostly considered (J, L)-regular QC codes 
that are described by a J x L polynomial parity-check matrix. 
However, one can construct many interesting ( J', L')-regular 
QC LDPC codes with a J x L polynomial parity-check matrix 
where J' ^ J and/or V ^ L. In this section we present a 
technique for constructing such codes, along with analyzing 
them. Of course, there are many other possibilities; we leave 
them open to further studies. 

Example 28. Let J' = 3 and V = 4. Consider the code C 
defined by the type-II (J', L')-regular polynomial parity-check 
matrix H'(x) G Fi, r> [x] JxL in © (with J = J' and L = L'), 
which for ease of reference is repeated here 



H'(x) 4 



Now, consider the code C" defined by the following type- 
I ( J', L') -regular polynomial parity-check matrix H"(x) £ 
F< r> [x] JxL (with J = 2 J' and L = 2L') 



X + x 2 





x 4 


X 8 


X 5 




x 10 + x 20 








x 25 +x 19 





X 7 + X 14 



H"(x) = 



V 


X 2 








j a 





x 8 





x 2 


X 1 











x 4 





x s 


X 5 





x e 







x 20 











X 5 





Q 

X 


x 20 


x w 














x 25 


x 19 








x~ 


x 14 








x 19 


x 25 








x 14 


x 7 



Here, in the transition from H'(x) to H"(x), 1 x 1-submatrices 
are replaced by 2 x 2-submatrices in the following way: 

"0 


A submatrix like [x a ] is replaced by the submatrix 

'0 



The submatrix [0] is replaced by the submatrix 







(or the submatrix 



A submatrix like [x a 



x 
x° 



(or by the submatrix 



x" 

x a 

is replaced by the submatrix 

x b x a 
x a x b 



). 



It can easily be verified that the Tanner graph of C" is a 
double cover of the Tanner graph of C. This easily follows 
from the fact that in each of the above replacements, the 
2 x 2-submatrices are incidence matrices of double covers of 
the graphs described by the 1 x 1-submatrices; see also the 
discussion in Appendix Q] Similarly, the proto-graph of C" is 
a double cover of the proto-graph of C. Therefore, the Tanner 
graph of C" is not only an r-cover of its proto-graph but also 
a 2r-cover of the proto-graph of C, 

Applying the bounds (Q} and (f2]i to H"(x) we obtain, 
respectively, d min (C") 80 and d min (C") 108. We 
can improve these bounds with the help of Lemma [30] in 
Appendix [Q which gives d mhl (C") 2d min (C) = 2-32 = 64. 



(Similarly, we can use Lemma [30] to derive the lower bound 

dmin(C") 55 d min (C') = 32.) 

We remind the reader that for r = 46, the code C has 
parameters [184, 47, 32], see Example [T2l It turns out that for 
the same value of r the code C" has parameters [368, 93, 56], 
i.e., in the transition from C to C" we have have significantly 
increased the minimum Hamming distance from 32 to 56. (Of 
course, C" is twice as long as C' .) Note that for this value of 
r the minimum Hamming distance of C" is relatively close to 
the upper bound predicted by Lemma [30] in Appendix [I] and 
by increasing r to values larger than 46 we could expect an 
increase in the minimum Hamming distance. 

The code with r = 46 has potential to have better minimum 
Hamming distance and a Tanner graph with higher girth than 
the [184, 47, 32] type-II code C and hence, much better than 
the type-I code C that we started with in Example [TT] It would 
however not be optimal with respect to the two bounds (Q]i 
and (f2]i. This is compensated by the easy way in which it is 
obtained from the type-II code, and by the new upper bound 
of 64 that can be deduced from the cover structure, giving 
thus concrete information on the parameters of this code. 
We computed the parameters of this code for r = 46 and 
obtained that C" is a [368, 93, 56] code. Therefore its rate 
is roughly the same as the rate of the [184, 47, 32] type-II 
code C . The two codes have also the same (3, 4)-regularity. 
However, the minimum Hamming distance is significantly 
increased from 32 to 56. We also note that by increasing r 
to values larger than 46 we could expect an increase in the 
minimum Hamming distance, (not more than the bound 64 set 
by the cover structure). 

The sum-product algorithm decoding performance of code 
C" is shown in Figure[3]and compared to a randomly generated 
(four-cycle free) (3, 4)-regular LDPC code of the same length 
and nearly the same rate. □ 

This process can now be repeated. Start with the type-I 
polynomial parity-check matrix H"(x) and form a type-II or 
type-Ill polynomial parity-check matrix H"'(x) from which 
a 2-cover type-I matrix can be derived as above. Another 
method is to start with a somewhat dense matrix that has a 
large minimum Hamming distance and to take sequences of 
2-covers using the above method until a type-I polynomial 
parity-check matrix is obtained or a large minimum Hamming 
distance. We leave this for further studies. 

This construction is just one of many possible ways to 
construct and analyze a (J', L') -regular QC LDPC code with a 
polynomial parity-check matrix that has J rows and L columns 
with J' ^ J and V ^ L. We leave it for future research to 
construct and analyze such codes. 

VIII. Conclusions 

We have presented two minimum Hamming distance upper 
bounds for QC LDPC codes, one based on the polynomial 
parity-check matrix and one on the weight matrix. Afterwards, 
we have investigated how these upper bound can be strength- 
ened based on the knowledge of Tanner graph parameters like 
girth. We have also constructed several classes of codes that 
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Fig. 3. Decoding performance of the [368, 93, 56] QC-LDPC code vs. a randomly generated (four-cycle free) [368, 92] LDPC code under sum-product 
algorithm decoding when transmitting over a binary-input AWGN channel. (For more details, see Example [28]) We observe the onset of an error floor of 
the word error rate at about 4.5 dB for the randomly generated (four-cycle free) LDPC code. A similar observation was made for other randomly generated 
LDPC codes with the same parameters. (Not shown in the plot.) 



achieve (or come close) these minimum Hamming distance 
upper bounds. 

In future work it would be interesting to establish similar 
bounds for the minimum pseudo-weight (for different chan- 
nels) of QC LDPC codes. (Some initial investigations in that 
direction are presented in [33].) 



Appendix A 
Proof of Lemma|6] 

Let S = {io, ii, ■ ■ ■ , ij} be the chosen size-( J+l) subset. 
In order to verify that c(x) is a codeword in C, we need to 
show that the syndrome s T (x) = H(x) • c T (a;) (in F^[a;]) is 
the all-zero vector. For any j G [J], we can express the j-th 
component of s(x) as follows 

s j( x ) = hj,i(x)ci(x) = hj,j(x) perm (H s \j(x)) 

ie[L] ieS 

= h iA x ) det (H s \i(a;)) , 
ies 

where in the last step we used the fact that for rings with char- 
acteristic 2 the permanent equals the determinant. Observing 
that Sj (x) is the co-factor expansion of the determinant of the 



|<S| x |5 1 -matrix 





hj,h 0) 


h 3,rj ( x ) 


ho,i (x) 




ho,ij{x) 


hi,i (x) 


h ldl( x ) 




hj-i,i (x) 


hj-i^ix) ■ 


■ hj-i tij (x) 



and noting that this latter matrix is singular (because at least 
two rows are equal), we obtain the result that s(x) = and 
that c(x) is indeed a codeword in C, as promised. 

Appendix B 
Proof of Theorem[8] 

We use the same notation as in Lemma [6] and Theorem [7] 
(and their proofs). Then, because the permanent is a certain 
sum of certain products, we can use the triangle and product 
inequality of the weight function to obtain 

w H (c(x)) = ]Pwt fperm(H s ^(o;))^ 
^perm((wt(H(x))) sv .) 
= Perm (A sv ) , (8) 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



13 



where in the last step we have used the definition A = 
wt (H(x)). 

With this, the result in (0 is nearly established, but not 
quite. Namely, it can happen that the codeword construction 
in Lemma [6] produces the all-zero vector c(x); in Eq. (0 
of Theorem [7] we properly took care of this case by using 
the min* operator instead of the min operator. However, 
unfortunately such an all-zero vector c(x) can yield a nonzero 
term X^ieS P erm (-A-S\i) m ® aR d we need to take care 
of this degeneracy. Our strategy will be to show that the 
minimum Hamming distance of the code is never larger than 
such a nonzero term and therefore, although this nonzero term 
appears in the min* operation in (0, it does not produce a 
wrong upper bound. 

Example 29. Before we continue, let us briefly discuss a 
polynomial parity-check matrix where the above-mentioned 
degeneracy happens, namely, let C be a code with polynomial 
parity-check matrix 



H(x)4 
and with weight matrix 



1 

x 

l + x 



1 



1 



/(*)' 
h(x) 



1111 wt(f(x)) 
1111 wt 

2 2 2 wt (h(x) 



where f(x), g{x), and h(x) are some arbitrary polynomials 
such that h(x) ^ f(x) + g(x). 

Let S — {0,1,2,3}. Clearly, the matrix Hs(:r) is rank- 
deficient because the last row of this matrix is the sum of 
the first two rows. This implies that all 3 x 3 sub-matrices 
of Hs(x) have determinant zero, i.e., that all 3 x 3 sub- 
matrices of Hs(x) have permanent zero, and so the codeword 
generating procedure in Lemma [6] yields the codeword c(x) = 
(0, 0, 0, 0, 0) for the above choice of the set S. However, the 
term in (0 is 



perm (Ag 



ies 



= ((2-1-1 + 2-1-1) + (2-1-1 + 2-1-1) + (2-1-1 + 2-1-1)) + 

((0-1-1 + 0-1-1) + (2-1-1 + 2-1-1) + (2-1-1 + 2-1-1)) + 

((0-1-1 + 0-1-1) + (2-1-1 + 2-1-1) + (2-1-1 + 2-1-1)) + 
((0-1-1 + 0-1-1) + (2-1-1 + 2-1-1) + (2-1-1 + 2-1-1)) 

= 36 . (9) 

If we can show that the minimum Hamming distance of the 
C is not larger than 36 then we are sure that the value of 
Xaes P erm (As\i) does not yield a wrong upper bound. We 
will do this by exhibiting a nonzero codeword c'(a;) with 
Hamming weight not larger than 36. 

In order to construct such a codeword c'{x), let H'(x) be 
the 2x5 sub-matrix of H(x) that consists of the first two 
rows of H(.t), and let A' be the 2x5 sub-matrix of A that 
consists of the first two rows of A. Because of the above- 
mentioned rank deficiency of 11$ (x), any vector c'(x) that is 
in the kernel of H'(a;) and that is zero at the fifth position must 



be a codeword in C. Now, applying the codeword generating 
procedure in Lemma|6]for H'(x) and for the set S' = {0, 1, 2} 
yields the codeword 



c'(x) 



0) 



(Because of the choice of S', it is clear that the fifth position 
of c'(a;) is zero. Moreover, and most importantly, because the 
matrix H^x) has full rank, c'(x) is a nonzero codeword^) 
This nonzero codeword yields the minimum Hamming dis- 
tance upper bound 

£ perm (A' s , v ) = (1-1 + 1-1) + (1-1 + 1-1) + (1-1 + 1-1) 

= 6. (10) 

Clearly, 6 is not larger than 36, and so the term 
J2ies P erm (As\i) implied by c(x) is of no concern. 

The fact that the value in ( TTOb is not larger than the value 
in (0 can also be seen more indirectly from the following 
observation. Namely, multiplying the expression in ( TlOb by 2 
(the weight of the element in the third row and fourth column 
of A) we obtain 

(2-1-1 + 2-1-1) + (2-1-1 + 2-1-1) + (2-1-1 + 2-1-1) , 

which is a sub-expression of (0. Because all summands in (0 
are positive, it is clear that the value in ( TTOb cannot be larger 
than the value in (0. □ 

In order to complete the proof of Theorem [8] let us 
generalize the observation that we have just made in the above 
example. Namely, let S be a subset of [L] with \S\ = J + 1 
and let c(x) be the codeword that is obtained by the codeword 
generating procedure in Lemma [6] for the set S. Before 
continuing, let us observe that c(x) is the all-zero codeword 
if and only if all J x J sub-matrices of Hs(:z:) have zero 
permanent, if and only if all J x J sub-matrices of H^x) 
have zero determinant, if and only if the matrix Hg(x) is 
rank-deficient. 

We want to show that the value of ^ ie5 perm^Asyj), 
if it is nonzero, is always an upper bound on the minimum 
Hamming distance of C. 

• Assume that Hs(x) has full rank. Then c(x) is a nonzero 
codeword and so ^2 ieS perm(A 5 \j) is a minimum Ham- 
ming distance upper bound because of the inequalities 
in ©. 

• Assume that Hs(x) has not full rank and that 
X; ie5 perm(A 5 \ i ) = 0. Then perm(A s \j) = for 
all i £ S. It follows that perm(H(x) l 5\ i ) = for all 
i G S and that c(x) is the all-zero codeword. Therefore, 
although c(x) is the all-zero codeword, this case is 
properly taken care of by the min* operator in (0. 

• Finally, assume that Hs(:e) has not full rank and that 
Eies P erm (-A-<S\i) > 0- (Note that this can only happen 
for J 2.) Without loss of generality, we can assume 
that the rows of H(.t) are ordered such that the last row 
of Hs(x) is linearly dependent on the first J — 1 rows 

7 In the proof of the general case we will also have to take into account the 
case where H^/fl) does not have full rank. 
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of Hs(x). Denote the entries of H(x) by hj^(x) and the 
entries of A by a^j. Moreover, let H'(x) be the (J — 
1) x L submatrix of H consisting of the first J — 1 rows 
of H(x) and let A' be the submatrix of A that consists 
of the first J — 1 rows of A. 

Because of the assumption Ylies P crm ( A s\i) > 0' there 
must be at least one i e 5 such that permYAsu) > 0. 
Using the co-factor expansion of the permanent of A.s\i> 
this implies that there is at least one i* € S\i such that 



aj_ M . ■ perm(A' f 



(S\i)\i*. 



> 



(11) 



Fix such an i* and let S' = S\i*. Assume for the mo- 
ment that H' s ,(x) has full rank. Applying the codeword 
generating procedure in Lemma [6] for the polynomial 
parity-check matrix H'(x) and the set S' we obtain a 
nonzero vector c'(x) which is in the kernel of H'(x). 
Because of the rank deficiency of H^fi) and because 
= for i' G [L] \ S' (and therefore c£, (x) = for 
i' G [L] \ S), the vector c'(x) must also be a codeword in 
C. Therefore, because c'(x) is a nonzero codeword, the 
minimum Hamming distance of C can be upper bounded 
as follows 



(C) < wn(c'(x)) < ^2 P crm ( A 5' 

i'eS' 



V 



< aj-i,i« • P erm ( A s'\i') 
i'eS' 



where we used the fact that the inequality 

in (fTTT i implies aj_i !t « 1. However, because 

■ Si'e5' P erm ( A S'\i') * s a su b- ex P ress ion of 
E lS 5 P crm (A <s \ 4 ) we obtain 



(c) < E pcrm ( A,s 



where we used the fact that J2ies P crm ( A SV') contains 
only non-negative summands. 

It remains to deal with the case where H^, (x) has not full 
rank. We can deal with this case with a similar procedure 
as above. Note that the above choice of i* ensures that 
there will be a suitable i'* G S'. 

Appendix C 
Proof of Corollary[9] 

We have to consider two cases. First, assume that the weight 
matrix A is such that there is at least one set 5' C [L] with 
\S'\ = J + 1 such that J2ies' P erm ( A S'\i) > 0- Using the 
fact that the weight matrix A of a type-I QC code contains 
only zeros and ones, we can conclude that for such an S' we 
have perm (A-s'\ i) ^ J* f° r a ll i G S' , which implies that 

perm(A s \ l ) < ^ P erm ( A S'\i) 
ies' 

( J + 1) • J! = {J +1)1 . 

ies' 

This is the upper bound that we set out to prove. 

Secondly, assume that the weight matrix A is such 
that for all sets S C [L] with \S\ = J + 1 it holds 



dmin(C) < min* 

SC[L] 

isi=j+i ies 



that X)iG5 P erm ( A 5\i) — 0- This, notably implies that 
perm (A^) = for all sets S and all i G S. (Parity-check 
matrices with such a weight matrix A are rather degenerate 
and in general uninteresting. However, we need to properly 
take care of this case too in order to verify the corollary 
statement holds for all possible type-I QC codes.) It follows 
that H(x) must be such that for all sets S C [L] with \S\ = 
J + 1 and for all i G <S it holds that perm (Hs\j(x)j = 0, 
i.e., dct (HsyK 21 )) = 0- This latter statement, however, is 
equivalent to the statement that H(x) does not have full 
row rank. The code C can therefore also be defined by a 
suitably chosen (J— 1) x i-submatrix of H(x). If J > 1 
then applying this corollary recursively to this (J— 1) x L- 
submatrix we obtain d m - m (C) ^ ((J— 1)+1)! = J!, which 
implies d m - m (C) ^ (J+l)!. Otherwise (i.e., when J = 1), it 
clearly holds that d min (C) < (J+l)! = 2! = 2. 

Appendix D 
Proof of Theorem[T71 

Assume that C has polynomial parity-check matrix H(x). 
We establish upper bounds on the girth of the Tanner graph 
of H(.t) by exhibiting the existence of certain cycles in that 
graph. These cycles will be found using techniques that were 
discussed in l29l , Q. 

a) Let 

r x a x x c 

d e f 

X X X 1 

be any submatrix of H(x) having the first weight 
configuration. The path 



x a — > x d — > — > x c — > x b — > x e — 
-> x c -> x s -> x e -> x b -> x a , 



x d -» x a 



shows that the Tanner graph of H(.t) has at least one 
12-cycle since 

(a-d) + (f-c) + (b-e) + (d-o) 
+ (c-/) + (e-6) = 

in Z (and therefore also in Z/rZ). 
b) Let 



be a submatrix matrix of H(x) having the second weight 
configuration. The path 



x a -> x b -> a; d 



.r 



shows that the Tanner graph of H(.t) has at least one 
10-cycle since 

(a-b) + (d-e) + (b-a) + (c-d) + (e-c) = 

in Z (and therefore also in Z/rZ). 
c) Let 

[x a + x b x c + x d ] 
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be a submatrix matrix of H having the third weight 
configuration. The path 



shows that the Tanner graph of H(x) has at least one 
8-cycle since 

(a-b) + (c-d) + (b-a) + (d-c) = 

in Z (and therefore also in Z/rZ). 
d) Let 

[x a + x b + x c ] 

be a submatrix matrix of H having the stated weight 
configuration. The path 

x a -> x b -> x c -> x a -> x b -> x c -> x a 

shows that the Tanner graph of H(x) has at least one 
6-cycle since 

(0-6) + (c-o) + (b-c) = 
in Z (and therefore also in Z/rZ). 

Appendix E 
Proof of Lemma [T9l 

Note that a 2x2 submatrix B(x) must, up to row and 
column permutations, look 



either like 
like 



or like 
or like 



x x 
x d 
h 



X 





or like 
or like 



x 


x a 








for some a,b,c,d £ Z/rZ. 
In the first case, 



wt (perm (B(x))J < perm ( wt (B(x) 
holds if and only if 



wt (x a+d + x b+c ) < 2 



if and only if 
if and only if 
if and only if 



x a+d + x b+c 



(inF< r> [x]), 



a + d = b + c (in Z/ rZ) , 



c + d - b = (in Z/rZ) , 



which is equivalent to the existence of a 4-cycle in the Tanner 
graph, see the conditions in 1291 . J5). 

In the second and third case, wt(perm(B(x))J < 
perm(wt(B(x))) holds if and only if wt(x a+d ) < 1, i.e., if 
and only if 1 < 1. However, this is never the case. This agrees 
with the observation that such 2x2 submatrices cannot induce 
a four-cycle in the Tanner graph |29l , J5). 

In the fourth, fifth, and sixth case, wt(pcrm(B(x))) < 
perm(wt(B(x))) holds if and only if wt(0) < 0, i.e., if and 



only if < 0. However, this is never the case. This agrees with 
the observation that such 2x2 submatrices cannot induce a 
four-cycle in the Tanner graph (29j, 0. 

The proof is concluded by noting that a 4-cycle can only 
appear in a 2 x 2-submatrix of a type-I polynomial parity-check 
matrix. 

Appendix F 
Proof of TheoremI2T1 

The proof of this theorem is based on upper bounding 
the minimum Hamming distance upper bound in ((T), thereby 
taking advantage of the fact that H(x) is assumed to be of 
type I and to have a 4-cycle. For the ease of reference, (HJ is 
repeated here, i.e., 

dmin(C) < nun* ^ wt (perm (H 5 \ ?; (x))) . (12) 

si=j+i ies 

Without loss of generality, we can assume that the rows and 
columns of H(x) are labeled such that the present 4-cycle 
implies that 

h Q0 (x)h n (x) + h Q1 (x)h 10 (x) = (in F< r> [x]) . (13) 

We define the following sets. 

• Let S^s{ be the set of all sets S C [L] with |<S| = J+l 
that are not supersets of {0, 1}. 

• Let S^jo.i} be the set of all sets S C [L] with \S\ = J+l 
that are supersets of {0, 1}. 

Then ( fT2l can be rewritten to read 



E 

(14) 

The first argument of the min-operator in ( fT4l can be dealt 
with a reasoning that is akin to the reasoning in the proofs of 
Theorem|8](cf. AppendixlBt and Corollary |9](cf. Appendix P. 
This yields 



dmin(C) < min min* V] wt ( perm (H S \i(x)) ) , 

(perm (H S \j(x)) 



min* y wt 



(perm (H SVi (x))) < (J+l)! . (15) 



Therefore, let us focus on the second argument of the min- 
operator in ( fT4l . i.e., 



ies 



wt ( perm (H s \j(x)) 



(16) 



We consider two sub-cases. First, assume that the polynomial 
parity-check matrix H(x) is such that there is at least one set 

S' £ Sd{os} such that E ie 5' wt (p erm ( H s'\»( :z: ))) > °- 
Any upper bound on this sum will be a valid upper bound on 
the expression in ( TTcl l. For i £ {0, 1} we find that 



wt (pcrm(H 5 ,\ 4 (x))) < J! , 



(17) 



where we used the fact that H(x) is a type-I polynomial parity- 
check matrix. 
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For i G S' \ {0, 1}, however, we want to use a more refined 
analysis for which we define the following sets. 

• We define V' to be the set of all permutation mappings 
from [J] to S' \ {%}. 

■ We define V" C V' to be the set of all permutation 
mappings from [J] to S' \ {i} that map (0, 1) to (0, 1) 
or map (0, 1) to (1,0). 

• We define V'L ^ to be the set of all permutation map- 
pings from [J] \ {0, 1} to S' \ {0, 1, i}. 

With these definitions, we obtain 

perm (H 5 /\j(x)) 

= e n h oMo)( x ) 

= e e n^>«o) 

<?ev> je[J] aev\v" je[J] 

= (h 0Q (x)hi 1 (x)+h 01 (x)h w (x)) ■ ^ [ h j^U)( x ) 

"CP'^^y j£[J]\{o,i} 

+ e n h w)i x ) 

*£V'\V" j£[J] 



aev\v"je[J] 



[ h jMj) (x) (in¥ { 2 r) [x}) 



where in the last equality we have taken advantage of (13[ . 
Clearly, \V'\V"\ = \V\-\V"\ = J! - 2(J - 2)!, so that we 
can upper bound the weight of pcrm(H < 5/\ i (x)) as follows 



wt (perm(H <s , v (x))) ^ J! - 2(J-2)! 



(18) 



where we have again used the fact that H(ir) is a type-I 
polynomial parity-check matrix. Combining (fTTI i and d!8t . we 
obtain 



E 



wt^pcrm (H s «y(3; 

wt (perm (H S /\j(x)) 



E 

ie{o,i} 



+ E wt f perm (H S /\j(i 

i65'\{0,l} 

2J! + (J-1)(J! - 2(J-2)!) 
==: (J+l)! - 2(J-1)! . 



(19) 



It remains to deal with the second sub-case, namely 
where we assume that the polynomial parity-check matrix 
H(x) is such that for all sets S G <Sd{o.i} ^ holds that 
Eies w * (p crm (Hs\i(a;))j = 0. This, however, is equiv- 
alent to the assumption that for all sets S G ^{o.i} an d a ll 
i G S it holds that perm (H 5 \j(x)) = det (H S u(x)) = 0, 
which in turn is equivalent to the assumption that for all 
S G Sd{o.i} the sub-matrix H^fi) does not have full row 
rank. 

Pick any S' G S'Djo.i} an d let code C be the code defined 
by Hs'. Without loss of generality, we can assume that the 
rows of H(a;) are ordered such that the last row of Hs> (x) is 
linearly dependent on the first J — 1 rows of H5/ (1). Let C" 



be code that is defined by the (J — 1) x ( J + 1) submatrix 
H[j_i], S ' of H s >. Clearly, 

eUn(C) < daun(C) = d min (C") < ((J-l) + l) ! = J!, 

where the first step follows from the fact that any nonzero 
codeword in C induces a nonzero codeword in C, the second 
step follows from the equivalence C and C", and the third 
step follows from Corollary [9] Without loss of generality we 
can assume that J ^ 2 (otherwise a type-I polynomial parity- 
check matrix H(a;) cannot have a four-cycle), and so we can 
upper bound the previous result as follows 

dmi„(C) (J+l)! - 2(J-1)! . (20) 

Finally, combining ( fT3T > and (|19t . or ( TT3T > and ( |20t , we 
conclude that the theorem statement is indeed correct, inde- 
pendently of which sub-case applies. 

Appendix G 
Proof of Lemma[231 

We consider only the case where all the entries of B(x) are 
monomials, i.e., 



B(ar) 



x 
x e 



X J 
X i 



for some a,b,c,d,e, f,g,h,i G Z/rZ. (The discussion for 
matrices B(x) where some entries are the zero polynomial is 
analogous.) By expanding the permanent pcrm(B) of B(z) 
we obtain 



x a+e+i +x b+ f +9 - 



X 



c+d+hj_ x c+e+g j_^,f+h+a 



X J 



„d+b+i 



Therefore, 



wt 



perm (B(x)) ) < perm ( wt (B(x)) 



holds if and only if 



x a + e + i ^ a; b +/+9 . 



-.7.: 



c+d+h j^^c+e+g j^ n ,f +h+a j^ n ,d+b+i 



-X 



< 6 



For this to hold, there must be at at least two monomials that 
are the same (in [x]). Two different cases can happen. 

• Suppose that two monomials like x° +e+l and x b+ l +g are 
the same (in [x]). Then 

a + e + i = b + f + g (in Z/rZ) , 

i.e., 

a-g + i-f + e-b = (in Z/rZ) . 

According to conditions in ||29l , f5| , this is equivalent to 
the existence of a 6-cycle in the Tanner graph. 

• Suppose that two monomials like x a+e+l and x a+ f +h 
T{r) [x]). Then 

= a + f + h 



are the same (in F 
a + e + 
i.e., 

e — h 



(in Z/rZ) , 



- / = (in Z/rZ) . 

According to conditions in ||29l , 0, this is equivalent to 
the existence of a 4-cycle in the Tanner graph. 
The proof is concluded by noting that a 6-cycle can only 

appear in a 3 x 3-submatrix of a type-I polynomial parity-check 

matrix. 
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Appendix H 
Proof of Theorem|241 

The proof is very similar to the proof of Theorem [21] in 
Appendix [F] and so the we will only discuss the main steps of 
the argument. Here, without loss of generality, we can assume 
that the rows and columns of H(x) are labeled such that 

h 00 (x)h n (x)h22(x) + h 01 (x)h 12 {x)h 2 o{x) = . 

The following steps are necessary to adapt the proof of 
Theorem [21] in Appendix [F] to the present corollary. 

• ^{o 1} i s replaced by a similarly defined set S^^ 1 2 y 

• £b{o,i} is replaced by a similarly defined set Szi{o,i.2}- 
> The line of argument leading to ([l9T l is replaced by the 

observation that for any S' e S^fcu^} an d any i £ S' it 
holds that 

wt (pcrm(H s , v (x))) < J! - 2(J-3)! . 
Therefore, for any S' € Sb-fci^} we nave 
^2 wt (perm(H,s'\;(x))) 

lES' 

= wt^pcrm(H s , v (j;)) 

i£{0,l,2} 

^2 wt (perm (H SI y(x)) 

iS5'\{04,2} 

^3J! + (J-2)(J!-2(J-3)!) 
(J+l)! - 2(J-2)! . 

• The line of argument leading to (|20T > is replaced by the 
observation that 

cUn(C) < daun(C) = d min (C") < ((J-l) + l) ! = J! . 

Without loss of generality we can assume that J ^ 3 
(otherwise a type-I polynomial parity-check matrix H(x) 
cannot have a six-cycle), and so we can upper bound the 
previous result as follows 

dmUQ < (J+l)! - 2(J-2)! . 

Appendix I 
Proof of Theorem|251 

The proof of the first part of the theorem is very similar 
to the proof of Theorem [21] in Appendix [F] and the proof of 
Theorem [24] in Appendix [H] and is therefore omitted. 

Therefore, let us focus on the second part of the corollary. 
Because the products on the left-hand and right-hand side 
of (0 are assumed to be nonzero, for every j e 1Z there 
exists an integer Pj, a (j) G Z/rZ such that 

h j,<r(j)( x ) = xP3 -" <j) ■ 

Similarly, for every j £ R there exists an integer Pj. T (j) £ 
Z/rZ such that 



The condition (0 can then be rewritten to read 

Y[ = Y[ x P3Mf> On F< r> [x]) . 

jen jen 

Clearly, this holds if and only if 

Y p j>« = Y p ^u) (in z / rZ) • (21) 

In order to continue, we define it to be the permutation 
mapping n = <t _1 o t from 1Z to 1Z. Moreover, we let j' be 
some element of 1Z and we set j' t = tx* (j' ), t = 1, . . . , R—l. 
Then, the condition in ( f2TT > holds if and only if 



R-l 



R-l 



Y PX,*Ul) = Y p A,r(j' t ) ( in Z M) . 
t=a t=o 

which holds if and only if 

Y (PiXii) -Pji,rUl)) = ([n Z / rZ > • W 

1=0 

Because we assumed that it = er -1 o r is a cycle permutation 
of order R, the condition in (l22l is equivalent to Tanner's con- 
dition of the existence of a 2_R-cycle, see 1291 , 0. (Note that 

<r(j't+i) = «Hft)) = r{j[) and that a{j' ) = o{v{j' R _ x )) = 

T(f R -l).) 

Appendix J 
Graph Covers 

This appendix collects some results that are used in Sec- 
tion En] 

Let C be a QC code with polynomial parity-check matrix 

U(x) = [h j , i (x)]..G¥^[x] JxL . 

We define the decomposition 

H(x) = HW(i) + H( 2 >(x) (in ¥ 2 r) [x] JxL 
with matrices 

H«(x) = \h$(xj\ e¥ 2 r) [x} JxL 



3,1 



and 



H^(x) 



e¥ { 2 r) [x} JxL 



Based on this decomposition, we define a new code C with 
the polynomial parity-check matrix 



H(x) 



H«(x) nW(x) 
H< 2 )(x) H^(x) 



G F< r> [x] 2Jx2i 



Lemma 30. The minimum distances of C and C satisfy 

d,„i„(C) < d min (C) < 2 • d min (C). (23) 

Proof: Let us start by proving the first inequality in d23l . 

Let c(x) = (c' 1 ) (x), c^ 2 ^ (x)) be a codeword in C with 
Hamming weight wn(c(x)) = iuh (c^ 1 ' (x)) +bjh(c' 2 '(j;)) = 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



18 



d min (C). Let us show that c(x) = cW(i) + c^(x) G C. 
Indeed, because c G C, we have (in Fi, r> [x]) 

H « (x) ■ c« (x) T + H< a > (x) ■ (x) T = T , (24) 
H^ 2 ) (x) ■ c« {x) T + H« (x) ■ c< 2 > {x) T = T . (25) 

Adding these two equations we obtain (in [x]) 

(hW(x)+H( 2 )(^)) • (cW(x)+c( 2 '(.t)) T = 0, 

showing that c(x) G C. If c(x) ^ then 

^min (C) < wn(c(x)) < wh (c (1) (x) + c (2) (x)) 

= wn(c(x)) = d min (C) , 

thus proving the first inequality in ( l23l l. If c(x) = then 
(x) = (x) and so both (124-b and d25l l can be rewritten 
to read (in F^, r> [x]) 

(h«(x)+H( 2 >(x)) •c«0r) T = O, 

showing that c (1) (x) = c (2) (x) G C. With this, 

rfmin(C) < (c (1) (a;)) < 2 ■ WH (c (1) (^)) 

= wn(c(x)) = d min (C) , 

thus proving the first inequality in ( |23| >. 

We now prove the second inequality in ( |23] >. Let c(x) be a 
codeword in C with Hamming weight u>h(c(x)) = d m [ n (C) 
and define c(x) = (c(x),c(x)). Let us show that c(x) G C. 
Indeed, because c G C, we have (in F< r> [x]) 

H(x) ■ c(a;) T = . 
Therefore (in F^ [x]), 

(h (1 >(x)+H( 2 '(x))c(x) t = 0, 

(r) 

and so (in F2 [x]) 

H (1) (x) • c(x) T + H< 2) (x) ■ c(x) T = 0, 

H^(x) • c(x) T + HW(x) • c(x) T = . 

which are exactly the equations that c(x) must satisfy in order 
to be a codeword in C. Therefore, because c(x) ^ 0, 

dmin(C) ^ U>b(c(x)) = 2 ■ W H (c(x)) = 2 • d min (C) , 

thus proving the second inequality in (l23l >. ■ 

Assume that the matrices H^-'(x) and H' 2 ^(x) are such that 
when hji(x) and h^(x) are added to obtain hj^(x), then no 
terms cancel for any j and and i. In this case it can easily 
be seen that the Tanner graph of H(x) is a double cover of 
the Tanner graph of H(x). This means that Lemma [30l relates 
the minimum distance of a Tanner graph and the minimum 
distance of a certain double cover of that Tanner graph. 

There is another way to obtain the same double cover (up 
to relabeling of the coordinates). Namely, based on code C we 
define the code C with the polynomial parity-check matrix 



This time the matrix H'(x) is obtained from H(x) by replac- 
ing, for each j and each i, the lxl entry 



h hl {x) = hf}{x)+hf}{x) 



by the 2 x 2 entry 



hfUx) ^ i ) (x) / 



It can easily be checked that H'(x) equals H(x) up to 
reshuffling of rows and columns. Therefore, the Tanner graphs 
of H(x) and of H'(x) are isomorphic, showing that they define 
the same double cover of the Tanner graph of H(x) (up to 
relabeling of the bit and check nodes). 

Let us remark that similar 2-covers as in this appendix 
were also considered in IT341 and (35). Namely, for some 
(scalar) parity-check matrix H, these authors first constructed 
the trivial 2-cover 



H 



H 
H 



Then, in a second step they split H into H = H' + H" 



whereby the matrices H' and H' 
nonzero entries do not overlap, 
modified 2-cover as follows 

H' 



H' 



H' 



were chosen such that the 
Finally, they formulated a 



H" 
H' 



H'(x)GFr[c 



(r) r 12JX2L 



We conclude this appendix by remarking that similar results 
can be proved for general M-covers for any M ^ 1. 
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